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Abstract. Models of complex networks are generally defined as graph sto- 
chastic processes in which edges and vertices are added or deleted over time 
to simulate the evolution of networks. Here, we define a unifying framework 
— probabilistic inductive classes of graphs — for formalizing and studying 
evolution of complex networks. Our definition of probabilistic inductive class 
of graphs (PICG) extends the standard notion of inductive class of graphs 
(ICG) by imposing a probability space. A PICG is given by: (1) class B of 
initial graphs, the basis of PICG, (2) class 1Z of generating rules, each with 
distinguished left element to which the rule is applied to obtain the right el- 
ement, (3) probability distribution specifying how the initial graph is chosen 
from class B, (4) probability distribution specifying how the rules from class 
Ti are applied, and, finally, (5) probability distribution specifying how the left 
elements for every rule in class 1i are chosen. We point out that many of 
the existing models of growing networks can be cast as PICGs. We present 
how the well known model of growing networks — the preferential attachment 
model — can be studied as PICG. As an illustration we present results regard- 
ing the size, order, and degree sequence for PICG models of connected and 
2-connected graphs. 



1. Introduction 

Recent surge of empirical results pertaining to the characteristics of complex 
networks, from technological networks (the Internet, the power grid) to networks 
of social contacts and protein interactions, has prompted theoretical research in 
models of complex networks. Existing theoretical results have identified two mech- 
anisms that result in synthetic networks with characteristics closely matching those 
of the real- world counterparts. Models of complex networks are usually defined as 
graph stochastic processes in which edges and vertices are added or deleted over 
time to simulate the evolution of networks, based on (1) preferential attachment 
mechanisms [Tj, (2) copying mechanism [3], or a combination of the two [TIT]. The 
result of a graph stochastic process is a class of random graphs, which can also be 
obtained via probabilistic inductive classes of graphs (PICG). 

2. Definition of PICG 

Eberhard [7] was the first to inductively define a class of graphs. A comprehensive 
survey of existing inductive definitions for classes of graphs is presented in [2J. 



Date: December 27, 2006 (submitted). 

Key words and phrases, inductive classes of graphs, probabilistic methods, growing networks. 

This work was completed with the support of the Department of Applied Mathematics, Faculty 
of Mathematics and Physics, Charles University, Prague grant MSM0021 620838; the DELIS-IST 
001907 project of the European Comission; by the Ministry of Higher Education, Science and 
Technology of Slovenia, Project Jl-6062; and partially by the Slovenian Research Agency, Project 
P5-0168. 



1 



NATASA KEJZAR 1 , ZORAN NIKOLOSKI 2 , AND VLADIMIR BATAGELJ 3 

R2:£ 

*- | 

R 




Example construction: 
Rl Q Rl 



Figure 1. An example of an ICG class X = (B;R1,R2) and a 
possible construction of a graph in that class. 

Definition 2.1. An inductive class of graphs (ICG) X = (B;1Z) is defined by [5]: 

(1) class B of initial graphs, the basis of ICG, 

(2) class TZ of generating rules, each with distinguished left element (part of a 
graph) to which the rule is applied to obtain the right element. 

The rules describe a substitution of a part of a graph (left side) by another part of 
a graph (right side) and result into a graph again. A class X consists exactly of the 
graphs that can be constructed from the graphs in the basis B by applying a finite 
number of generating rules TZ. Figure [T] shows an example of an ICG X = (B;TZ), 
TZ = {Rl, R2} and a possible construction for one of its graphs. 

The small triangles attached to the nodes in the pictorial depiction of the rules 
presented in Fig. [T] denote or more edges. 

Here we extend the classical notion of inductive class of graphs by introducing 
a probabilistic space in the following manner: 

Definition 2.2. A probabilistic inductive class of graphs (PICG), X, is given by: 

(1) class B of initial graphs, the basis of PICG, 

(2) class TZ of generating rules, each with distinguished left element to which 
the rule is applied to obtain the right element, 

(3) probability distribution specifying how the initial graph is chosen from class 
B, 

(4) probability distribution specifying how the rules from class TZ are applied, 
and, finally, 

(5) probability distribution specifying how the left elements for every rule in 
class 1Z are chosen. 

The probabilistic inductive class X is closed under the class of generating rules 
1Z, i.e. a random graph G is in the PICG X if and only if it can be obtained by 
applying the generating rules in the given class 1Z finite number of times, starting 
with an initial graph from the basis B. All steps in this construction have to have 
positive probabilities. 

It is often desired that the generating rules are simple enough not only to facil- 
itate rigorous analysis but also to capture the mechanism for generating a class of 
graphs that match the properties of real- world networks. Two such properties are: 

(1) locality: the rule is local if and only if its left element is connected. 

(2) expansion: the rule is expanding if a selected property of graphs is increased 
(or remains the same) by the application of the rule (for example number 
of vertices, number of edges, girth — the size of the shortest cycle, etc.) 
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In general it can happen that the selected rule can not be applied in a current 
graph. In this paper we shall limit our discussion to the simple PICG definitions 
in which the left element of the selected rule can always be found. 

In order to prove a certain result for ICGs (as well as PICGs), the principle of 
inductive generalization is often used: 

Principle: Inductive generalization 

Given an ICG T and a property of graphs Q. If it holds that: 

(1) each graph in B has the property Q, and 

(2) each generating rule in 1Z preserves the property Q — if the graph on which 
the rule is applied, has the property Q, then also the resulting graph has 
the property Q, 

then all graphs from X have the property Q. 

An example: If all the elements (initial graphs Bpicg arid the rules IZpicg) 
in the definition of PICG have positive probabilities, then Xpicc — ^-iCGi where 
tiCG = (&PICG', T^-picg), because every graph from ICG has a positive probability 
to be created. 

Notation: Given a graph G = (V, E), let us denote the number of vertices 
|V(G)| = n, and the number of edges |2?(G)| = m; and random variables for the 
number of vertices with n and m for the number of edges. Let d(u) denote the 
degree of vertex u G V (G) — number of edges with u as an endpoint. 



3. An illustration: Preferential attachment model as PICG 

The model was presented by Barabasi and Albert in 1999 pQ. Its idea is to 
capture the evolution of the network that can possibly explain the emergence of a 
power-law degree distribution. The model is based on the following concepts: (1) 
a network grows in time (the number of vertices and number of edges grow) and 
that (2) a newcomming vertex is more likely to build an edge to a more "popular" 
vertex in the network (a vertex with relatively large number of edges). Barabasi 
and Albert called this principle a preferential attachment (abbreviated as pa). 

The algorithm of a model consists of two processes in a time step: an addition of a 
new vertex and a creation of m pa edges from this vertex to already existing vertices, 
which are chosen proportional to their current degree. To make the formulation 
of the model precise, consider to start the process with two vertices linked by 
m pa paralel edges. When adding a new vertex, the edges will be added one at a 
time (so for the second and subsequent edges the probabilities will be calculated 
with updated vertex degrees). This formulation has a nice property that only the 
m pa — 1 model needs to be analyzed [3]. To get the models with m pa > 1, vertices, 
added at time steps jm pa , jm pa — 1, • ■ ■ , (j — l)m pa + 1, should be collapsed into a 
single vertex j (added at time step j). 

We can describe this model as a probabilistic ICG very easily: The basis consists 
of an edge linking 2 vertices. The rule is just one — to a selected vertex an edge 
that ends with a new vertex is added (see Fig. [2|). We only have to specify, how 
the left element (the vertex) is chosen. The probability of selecting a vertex, where 
to apply the rule, is proportional to its current degree. The probability of a vertex 
i to be chosen is exactly Si = d(i)/Y^,j d(j). 
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Figure 2. The basis, the only rule for a preferential attachment 
model and the graph after the first step (after one application of 
the rule). 



4. Order and size of PICG 

For the simple PICG it is easy to determine the expected order (number of 
vertices) and expected size (number of edges) of the graph after t steps. We divide 
the rules in three different classes: expanding, shrinking and stable class. The 
change in the number of elements (vertices or edges) after the application of the 
rule in these classes is either growing, decreasing or staying the same. In order 
to get non-trivial results for large t, the rules' expected change in the number of 
vertices or edges has to be positive. Otherwise the graph dies out when t — > oo. 

Given a graph G G X, let p t (n = n) denote the probability that at step t the 
order of G is n. Consider Ri G 1Z a fixed chosen rule in X, Arij the number of 
vertices, that rule Ri adds to the graph and rj the probability of selecting rule Ri. 
The probability p t (n = n) can be expressed recursively as 

p t (n = n) = ^2 nPt-i (» = n - Am), 

with initial values for all basic graphs Bj, G B: po (n = k) = J2Bi-\Bi\=k Qi> where Qi 
is the probability, that Bi is chosen for initial graph. For the clarity of the further 
equations the following abreviation will be used p t (n) =pt(n = n). 
The expected number of vertices is therefore 

oc 

E t [n]=J2m (i) 

8=0 

The expected change in the number of vertices is then easily computed as 
(4.1) Hm ^1 = ]T nAn, 

Proof. We can write 

oo 

S *[ n ] = X/ r jPt -i (i - Anj) 

i=0 RjeK 

oo 

= H rj5^ipt_i(»-Anj) 

Rj en i=o 
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We split the inner sum into two sumands 

oo oo oo 

ipt-i (i - Anj) = (i - Anj)p t -i (i - An 3 ) + ^ ArijPt-i (i - Anj) 

i=0 i=0 i=0 

oo 

= kpt-i (k) + Arij 

fc=0 

= E t _i[ii]+Arij. 
Now, the whole equation can be written as a recursive relation 

E t [n] = »V [n] + An,-) 

= Ek-i [n] + rjArij. 

with a solution 

E t [n] = (t — k) r j An j + E k [n] 
from which the result follows. 

□ 

The expected number of edges can be calculated similarly, one just has to change 
the elements that one counts. Therefore 

oo 

E t H = X/ Pt ( m = 

j=0 

and the expected change in the number of edges can be calculated as 

lim — j" 1 ^ = riArrii. 

Rien 

The calculation of expected number of vertices and edges in a step t for a prefer- 
ential attachment model (see section [3]) is straightforward. Because of the applica- 
tion of a single rule, we know the exact numbers (the change of vertices and edges 
in a single time step is 1), and since after the first step t = 1, we have already 3 
vertices and 2 edges, we obtain 

E t [n pa }=t + 2 and E t [m pa ]=t + 1. 

In the next three sections, we present three PICG models. The ICG models are 
defined by (1) subset of the basis graphs, and (2) the subset of rules from Fig. [3] 
Rules are Rl — add a vertex with an edge to an arbitrary vertex with degree i > 0; 
R2 — add an edge to two vertices with arbitrary degrees (the vertices of the left 
side are different, but can be endpoints of the same edge, thus the rule allows for 
multiple edges but no loops); R3 — expand an edge with a vertex and another 
edge; R4 — add two interconnected vertices to a vertex with arbitrary degree. 
They share some common characteristics. (3) The basis of all models considered 
in this paper contains only one graph, therefore it is chosen with probability 1; (4) 
each of the rules in a model is chosen with a non-zero probability and; (5) the left 
elements from the rule are chosen uniformly at random. 

We then show how one can derive some probabilistic properties of graphs that 
belongs to the defined classes. 
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Figure 3. Bases (Si, Si?) and rules for the PICG models presented. 

5. Connected undirected graphs 

The ICG description is Z(S1; Rl, R2) according to the rules and the basis from 
the Fig. [3] For PICG extension, Tq, we add the probability for rule selection: Rl 
is selected with probability q and R2 with probability r = 1 — q, q <E (0, 1). 

5.1. Order and size. In the very precise formulation of the model we note, that 
in the first step only Rl can be applied, but the initial condition does not depend 
on the further evolution of the graph. 

The number of edges at step t is m if and only if the number of edges at step 
(t — 1) is (m — 1) and any rule is applied, therefore: 

Pt (m = m) = qpt-i (m = m — 1) + rpt-i (m = m — 1) = 
= pt-i (m = m - 1) . 

Lemma 5.1. The probability that a graph G € Ic has m edges at step t is given 
by: 



, , f 1 t = m 

Pt (m = m) = < , . 
y v ; | t ^ m 



Proof. The lemma follows by applying the initial condition, po (m = 0) = 1, to the 
previously stated recurrence. □ 

The number of vertices at step t is n if and only if the number of vertices at step 
(t — 1) is: (1) n and rule Rl is applied, (2) (n — 1) and rule R2 is applied, or (3) 
n and neither rule SI nor rule S2 is applied (thus, the probability is 0). 

Lemma 5.2. The probability that a graph G 6 Tc has n nodes at step t is given 
by: 

, \ f (* t ~ii)« n ~V- n+1 n>l 
Pt(n ) = < ^-n+l)^ . 
l_ n = or 1 

Proof. For the proof we have that: 



Pt (n) = qpt-i {n - 1) + rp t _i (n) . 

Let P n (x) = ^2 t Pt (n) . Multiplying the last expression by x and summing 
over all t, we get: 
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P n (x) = xrP n (x) + xqP n - X (x) , 

or equivalently, 

Pn (X) - j^—Pn-l (X) . 

1 — xr 

To determine the initial condition, i.e. P 2 (x), we find that: 

p t {2) = rp t . 1 {2)=r t -\ 
since p\ (2) = 1 and for every t > 0, p t (n = n) = when n = or 1. Therefore, 
-F2 (z) = i- rx - To obtain the desired result, we observe that p t (n) is the coefficient 
of the x l term in the expansion of P n (x): 



Pt (n) = [x^Pni 

= q n - 2 [a;*-™+ 1 ] 



n-2 



1 — xr 

1 



(l-xr) ( "- 2)+1 
2 , (t-„+l) + („-2)\ ,_„+, 



; t-n+1 



t-l 
t-n+1 



g n-2 r t-n+l i 



and we have the claim. 



□ 



Finally, we determine the expected number of nodes at step t. The result is 
summarized in the following theorem and corollary: 

Theorem 5.3. Given a graph G G Ic, the expected number of nodes at step t, 
denoted by E t [n], is: 

Et[n] = (t-l)q + 2-(t+l)q t - 1 . 

Proof. The proof involves some technical manipulation of the expression for the 
expectation: 



* ' t-l 



,»-2_t-»+l _ 



- E(i+ 2 )(,!;' 1 )^ t - j - i = 

= Di + 2 )(t''li) ( f irt " H -(' + 1 )« 



3=0 

t-l 



j=0 



Ej( t ! ? -li)^" i_1 + 2 -(*+ 1 )« 



t-j 



s 
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Since X)i=o-? (t-i)? Jr * 3 = ^ (^ ne expected value of the binomial distribution), 
we have the proof. 

□ 

Corollary 5.4. When t — » oo, Et ^ — > g m probability. 

Proof. Finding the limit gt j n l when i — * oo is an elementary exercise. □ 

5.2. Degree distribution. Degree distribution of several growing graphs can be 
approximately calculated with the use of generating functions [12] . Consider pd the 
probability of a vertex having a degree d. If there are n vertices currently in the 
graph, there are on average npd vertices of degree d currently in the graph. With 
this we assume, that the number of vertices of degree d is tightly concentrated 
around its expectation (npd), which corresponds with the simulation results for 
large n (see Fig. |4}. A claim similar to this one can be found rigourosly proven in 
the paper of Cooper et al [4] for scale free random graph processes. 

In the next step, the probability pd is going to be updated regarding the use of a 
rule. We denote it as p' d . The number of all vertices in the next step will be larger 
for the probability, that a new vertex is added, so n + q (only Rl adds a vertex). 
And the number of vertices with degree d present in the next step will be (n + q)p' d . 

In this way, we can write the rate equation for the evolution of the graph degree 
distribution. The number of vertices with degree d changes if either of the rules 
is applied. With probability q, the Rl is used, a vertex with degree 1 is added 
(linked to a randomly selected vertex). In the equation we describe the addition of 
this vertex by Kronecker Sdi ■ The choice of a vertex can influence pd iff the chosen 
vertex is of degree d (pd is diminished) or of degree (d — 1) (pd increases). 

R2 is applied with probability r. The vertices for rule R2 are chosen uniformly 
at random. When a vertex with degree (d — 1) is chosen, this adds to the number 
of vertices with degree d, but when a vertex with degree d is chosen, the number 
of vertices with degree d will be diminished. The probability, that two vertices of 
degree d are selected is 

np d (np d -\) , 1 \ 
2 _ Pd(npd ~ 1) _ 2 



»(»-!) n — 1 

2 



Pd- 



The entire change in the number of vertices with degree d in one time step can 
therefore be written as 

2 Pd-i _ 2 Pd +Pd-i(l -Pd-i - Pd) +Pd{l -Pd-i - Pd) =Pd-i -pd + Pd-i - Pd- 

The contributions, where both vertices change the pd probability simultaneously 
vanish in the limit of large n and have therefore been neglected. The rate equation 
follows: 

(n + q)p' d = np d + q (5 dl + p d -i - Pd) + r (pd-i - Pd) ■ 
p' d is taken to equal p d , when n — > oo, which gives 

(1 + q)p d - qSdi - Pd-l = 0. 

The solution of this equation can be found in terms of generating functions. Mul- 
tiplying by z d and summing over d (considering ^2 d Pd zd = G{z)) 

(1 + q)G(z) - qz - zG(z) = 

G ( z ) = 7T^\ ■ 

(1 + 9) -z 
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The [z'j-th term of the generating function G{z) corresponds to the probability pi 
and rewriting the generating function in terms of geometric series we get 

GW(z), _ q 



W) = 



lz=0 



(i + q) 1 



This shows, that the degree distribution of the graph follows the exponential func- 
tion. 

Considering the rule R2 only, where vertices in the left element are not linked, 
the inductive description X(B1;RI,R2) represents exactly the class of all simple 
connected undirected graphs. Note, that in this case, the left elements for R2 are 
not selected uniformly at random anymore. However this anomaly disappears, when 
t — » oo, while the probability, that the two random selected vertices share an edge 
goes to 0. 

The model above is therefore a good approximation for the class of all simple 
connected undirected graphs. This can also be seen from the simulation of the 
degree distribution (see Fig. |4|). 



i i i i i i i i i i i i i i i i 

1 3 5 7 91011 1415 17 



20 23 2526 2930 32 35 



degree 



FIGURE 4. 10 realizations of the model were grown for q = r = 0.5 
till they reached 10,000 vertices. The degree density was calcu- 
lated. Boxplot for each degree is represented on the graph (in 
black). We can see that variance of the degree densities is very 
small. Gray points represent the the mean-field degree distribu- 
tion fit. The lines are there just for easier visual inspection. 



6. 2- VERTEX-CONNECTED GRAPHS 

A graph G is 2-connected if \G\ > 2 and G — u is connected for every u £ V(G). 
In other words, no pair of vertices can be separated by the removal of any other 
vertex [BJ. The ICG description is T(B2;R2,R3) according to the rules and the 
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bases defined above [IT]. In the probabilistic extension Tiy, R2 is selected with 
probability q G (0, 1), and R3 with probability r, where q + r = 1. 

6.1. Order and size. The calculations for order and size are done very similarly 
to the ones in subsection 15.11 The change in the number of vertices and in the 
number of edges with the application of the rules is the same. The number of edges 
changes for 1 in every step and the number of vertices does not change when R2 
is applied (with probability q) and it changes for one when R3 is applied (with 
probability r). Note that the probabilities of choosing a rule are opposite to the 
ones in the subsection 15.11 The second difference comes with the different initial 
graph; B2 has already 3 vertices and 3 edges. 

Lemma 6.1. The probability that a graph G G X^y has m edges at step t is given 
by: 



pt (m = m) 



1 t = m + 3 
+ 3 



Proof. The lemma follows by applying the initial condition, po ( m = 3) = 1, to the 
recurrence stated in subsectionl5.11 □ 



Lemma 6.2. The probability that a graph G G Zijv has n nodes at step t is given 
by: 



Pi 



(n) J (t_,t + > n -y- n+3 n > 2 
[ ' \ ne {0,1, 



2} 



Proof. Similarly to the proof in subsection 15.11 



p t (n) = qpt-i (n) + rpt-i (n - 1) . 
For the generating function's expression, when P n (x) = ^2 t pt (n) x l we get: 

xr 

P n {x) = xqP n (x) + xrP n _ 1 (x) = P„_i (x) . 

1 — xq 

To determine the initial condition, i.e. P3 (x), we find that pt (3) = qpt-i (3) = 
g', since po (3) = 1 and for every t > 0, pt (n) = when n G {0, 1, 2}. Therefore, 
P 3 ( a; )= 1 . 

o \ / 1 — qx 

Again, pt (n) is the coefficient of the x l term in the expansion of P n (x): 

Pt{n) = ( 1 V-y- n + 3 , 

\t — n + 3/ 

and we have the claim. 

□ 

The expected number of nodes at step t is summarized in the following theorem: 

Theorem 6.3. Given a graph G G I2V > the expected number of nodes at step t, 
denoted by E t [n] , is: 

E t [n] = tr + 3 - ( * + ^ gV" 2 - (t + 2)tqr t - 1 - (t - 3)r*. 
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Proof. The proof involves some technical manipulation of the expression for the 
expectation: 



*n - g< i _, t +3 )»- i+ v-»=|:o + 3)( t !.),- 



3 r 3 



E (i + 3) ( * Vv - {t+l) f 1} g V- 2 - (t + 2)V- 1 - (t - 3)r* 
E i ^ * .)g*"V + 3 - {t + 1) f ~ 1} gV- 2 - (t + 2) V- 1 - (t - 3)r*. 



Since Y^j^oiit-j)^^ ^ = ^ r (^h e ex P ec ted value of the binomial distribution), 
we have the proof. 

□ 

Corollary 6.4. When t — > oo, 5^-1 — > r m probability. 

6.2. Degree distribution. We use the approximate, mean-field approach (as de- 
scribed in 15. 2p , solved by the generating functions methodology. The rate equation 
for the change in the number of vertices of degree d is written in the following way 

(n + r)p' d = np d + r (S d2 ) + q (p d -i - Pd) ■ 

p' d is taken to equal pd, when n — ► oo, which gives 

Pd - rd d 2 - qp d -i = 0. 

Multiplying by z d ~ 2 and summing over d (considering ^2 d >2Pd zd ~ 2 = F(z) — left 
shifted generating function, where po — and p\ = 0) 

F(z) - r - qzF(z) = 

T 

F(z) 



1 — qz 



The \z l 2 ]-th term (i > 2) of the generating function F(z) corresponds to the 
probability p% , 

J?(i- 2 )( z \ 

[z*- 2 }F(z) = -, m n = rq l ~ 2 . 

1 i y ' (i-2)\ lz=Q 

Figure [5] shows the fittings of degree distribution to the simulated data. 



7. 2-EDGE-CONNECTED GRAPHS 

G is called 2- edge- connected if \G\ > 1 and G — e is connected for every e £ E{G) 
[B]. The ICG description is T(B2; R2, R3, RA) according to the rules and the basis 
defined above [TT]. For the PICG model I2E, R2 is selected with probability q, R3 
with probability r and R4 with probability s, q + r + s = 1. 
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Figure 5. 10 realizations of the model were grown for q = r = 0.5 
till they reached 10, 000 vertices. Boxplot for the density of each 
degree is represented on the graph (in black). Gray points represent 
the the mean- field degree distribution fit. The lines are there for 
easier visual inspection. 



7.1. Order and size. The number of vertices at step t is n, n > 4, if and only if 
the number of vertices at step (t — 1) is: (1) n and rule R2 is applied, or (2) n — 1 
and rule R3 is applied, or (3) (n — 2) and rule R4 is applied: 

p t (n) = qp t -i(n) + rp t -i(n - 1) + sp t -i(n - 2). 

Consider P n (x) — X)t>oP*(") x *- Multiplying the last expression and summing 
over t > gives 

P n {x) = qxP n (x) + rxP n -i(x) + s.tP„_ 2 (x), 

or equivalently: 

(1 - qx)P n (x) = rxP n - 1 (x) + sxP n - 2 (x). 

Note that Pq(x) = P\{x) = Pi{x) = 0. From the initial condition, wc further 
have that Po(3) = 1, and pt(3) = q l po(3) = g*, therefore, Ps(x) = ■ 
The following claim can be inductively proven: 

Lemma 7.1. Given a positive integer n, the following holds for the generating 
function P n (x): 

• P Q (x) =Pi(x) = P 2 (x) =0, 

• M*) = Tk-v 

• For every n>3, P n {x) is a summation of [^] — 1 functions: 
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Pn(x)= ]T 
j=LfJ 



" 2 '2j - n + 2\ ^'-"+1^-3-2^-1 



n — j — 2 J (1 — 



Then we have the following theorem: 

Theorem 7.2. TTie probability that a 2-edge connected random graph has n vertices 
after applying t rules is given by: 



*(»>-e(1-;- + .%J + i)' 



3=LfJ 

Proof. Let c# = ( 2 ^"+ 2 2 )- Note that 
p t (n) = [x t ]P n {x) = 



t-j+l r 2j-n+l s n-3'-2 _ 



" 2 r 2i-n+l s n-i-2 a ,i-l 



n-2 

„n T ,2j-n+l n-j-2^t-i+l 1 1 



a^r J s J \x 



3-LfJ 



n-2 



\ " n 2j-n+l n-3-2/ * ) n t-j+l 



y < 



i=LfJ 



t-i+l r 2j-n+l s n-3-2_ 



□ 



According to the approximate calculation of the expected change in the number 
of vertices (Eq. (|4.1[) ). the expected change for 2-edge-connected graph model is 
— ^ sa r + 2s. 

The calculation of the number of edges involves even more demanding and te- 
dious calculations. For the clarity of the presented paper, we do not include that. 
According to Eq. (|4.1j) the expected change in the number of edges, when t — > oo, 
approximately equals to (q + r) + 3s. 

7.2. Degree distribution. The mean-field approach gives the rate equation for 
the change in the number of vertices of degree d 

(n + r + 2s)p' d = np d + r (S d2 ) + q (pd-i - Pd) + s (j>d-2 + 2S d 2 ~ Pd) ■ 

p' d is taken to equal pd, when n — > oo, which gives 

(1 + 2s)p d - (r + 2s)5 d 2 - qpd-i - spd-2 = 0. 

Multiplying by z d ~ 2 and summing over d (considering X)d>2P<2Z d ~ 2 = F( z )' a g am 
Po = and p\ = 0) 

(1 + 2s)F{z) - (r + 2s) - qzF(z) - sz 2 F(z) = 

— r+2s 
r \ Z > — (l+2s)-qz-sz i ■ 



14 



NATASA KEJZAR 1 , ZORAN NIKOLOSKI 2 , AND VLADIMIR BATAGELJ 3 



The [z l ~ 2 ]-th term of the generating function F(z) corresponds to the probability 
Pi. It can be expressed in a closed form according to Graham et al. [8] with the 
partial fraction expansion. 

We can write F(z) = P(z)/Q(z) and express Q(z) = (1 + 2s)(l — piz)(\ — p2z), 
because the coefficients of the generating functions of the form 1/(1 — az) are terms 
of simple geometric series expansion. Calculate the two roots of the "reflected" 
polynomial Q R (z) (pi,2), that represent the roots of the above expression of Q(z). 
[z n ]-th term of F(z) can then be expressed as a\p\ + a^p^ (when the roots are 
distinct). Oj is expressed as —pi(P(l/pi))/(Q'(l/pi)). 

In our case, the roots of Q R {z) equal 



_ q ± yV + 4s(l + 2s) 

Pl - 2 ~ 2(1 + 2s) 

The [z J ]-th coefficient of F(z) therefore is 



[z }F(z) = ai • v , K-s + a 2 - 1 



where 



2(1 + 2s) / \ 2(1 + 2s) 



(q+ vV + 4s(l + 2s))(r + 2s) 

a\ = , — ana 

2(l + 2s) V V 2 + 4s(l + 2a) 



a 2 



(q - ^/q 2 +4s(l + 2s)){r + 2s) 



2(1 + 2s)vV + 4s(l + 2s) 
Figure [S] shows the fittings of degree distribution to the simulated data. 

8. Conclusion 

In the paper we proposed a quite broad framework for the description of sta- 
tistically evolving complex networks, the probabilistic inductive classes of graphs. 
We presented, how the already existing network models can be described in this 
way. We give a general recursive equation for order and size of the classes of graphs, 
which belong to the simple PICGs and apply them on three PICG models of growing 
networks. 
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